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Abstract 

In this article, we deal with a multiple dimensional coupled Markovian BSDEs 
system with stochastic linear growth generators with respect to volatility 
processes. An existence result is provided by using approximation techniques. 
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1 Introduction 

Backward stochastic differential equations (BSDEs) was proposed firstly by Bismut 
(1973) [3] in linear case to solve the optimal control problems. Later this notion 
was generalized by Pardoux and Peng (1990) [21] into the general nonlinear form 
and the existence and uniqueness results were proved under the classical Lipschitz 
condition. A class of BSDE is also introduced by Duffie and Epstein (1992) [8] in 
point of view of recursive utility in economics. During the past twenty years, BSDEs 
theory attracts many researchers’ interests and has been fully developed into various 
directions. Among the abundant literature, we refer readers the florilegium book 
edited by El-Karoui and Mazliark (1997) [9] for the early works before 1996. Surveys 
on BSDEs theory also includes [10] which is written by El-Karoui, Hamadene and 
Matoussi collected in book (2009) [7] (see Chapter 8) and the book by Yong and 
Zhou (1999) [23] (see Chapter 7). Some applications on optimization problems can 
be found in [10]. About Other applications such as in field of economics, we refer 
to El-Kaoui, Peng and Quenez (1997) [11] . Recently, a complete review on BSDEs 
theory as well as some new results on nonlinear expectation are introduced in a 
survey paper by Peng (2011) [22]. 

One possible extension to the pioneer work of [21] is to relax as much as possible 
the uniform Lipschitz condition on the coefficient. A weaker hypothesis is presented 
by Mao (1995) [19] which we translate as follows: for all y,y,z,z and t E [0,T], 
the generator of the BSDE satisfies \f(t, y, z) — f(t , y, z) | 2 < n\y — y | 2 + c\y — y | 2 
a.s. where c > 0 and k is a concave non-decreasing function from R + to R + such 
that «(0) = 0, k(u) > 0 for u > 0 and f Q+ 1 /n(u)du = oo. An existence and 
uniqueness result is proved under such condition in [19]. Hamadene introduced in 
(1996) [13] a one-dimensional BSDE with local Lipschitz generator. Later Lepeltier 
and San Martin (1997) [18] provided an existence result of minimum solution for 
one dimensional BSDE where the generator function / is continuous and of linear 
growth in terms of (y, z). When / is uniformly continuous in z with respect to (to, t) 






Mu and Wu 


Page 2 of 14 


and independent of y , a uniqueness result was obtained by Jia [17]. BSDEs with 
polynomial growth generator is studied by Briand in [4] . The case of 1-dimensional 
BSDEs with coefficient which is monotonic in y and non-Lipschitz on 2 is shown in 
work [5]. About the BSDE with continuous and quadratic growth driver, a classical 
research should be the one by Kobylanski [12] (2000) which investigated a one¬ 
dimensional BSDE with driver | f(t,y,z) < (7(1 + |y| + | 2 | 2 ) and bounded terminal 
value. This result was generated by Briand and Hu into the unbounded terminal 
value case in [6] (2006). 

There are plenty works on one-dimensional BSDE. However, limited results have 
been obtained about the multi-dimensional case. We refer Hamadene, Lepeltier and 
Peng [14] for an existence result on BSDEs system of Markovian case where the 
driver is of linear growth on (y, 2 ) and of polynomial growth on the state process. 
See Bahlali [1] [2] for high-dimension BSDE with local Lipschitz coefficient. 

In the present article, we consider a high dimensional BSDE under Markovian 
framework as follows: 



for i = 1,2,...,n, with process X as a solution of a stochastic differential equa¬ 
tion (SDE for short). For each i = 1,2the coefficient Hi is continuous on 
(y 1 ,..., y n , z 1 ,..., z n ) and satisfies: 



As the generator of the BSDE, actually, Hi is of stochastic linear growth on Z % , or in 
another word, it is of linear growth w by to. Similar situation was considered in [15] 
in the background of nonzero-sunr stochastic differential game problem. However, 
in [15], the generator Hi is independent on (y l ,..., y n ). According to our knowledge, 


this general form of high dimensional coupled BSDEs system with stochastic linear 


growth generator has not been considered in literature. This is the main motivation 
of the present work. 

The rest of this article is organized as follows: in Section 2, we give some nota¬ 
tions ans assumptions on the coefficient. The properties of the forward SDE are also 
provided. The main existence result of BSDEs is proved in Section 3 where a mea¬ 
sure domination result plays an important role. This domination result holds true 
when we assume that the diffusion process of the SDE satisfies the uniform elliptic 
condition. For the proof of the main result, we adopt an approximation scheme fol¬ 
lowing the well know mollify technique. The irregular coefficients are approximated 
by a sequence of Lipschiz functions. Then, we obtain the uniform estimates of the 
sequence of solutions as well as the convergence result in some appropriate spaces. 
Finally, we verify that the limit of the solutions is exactly the solution to the original 
BSDE which completes the proof. 


2 Notations and assumptions 

In this section, we will give some basic notations, the preliminary assumptions 
throughout this paper, as well as some useful results. Let (S2, X, P) be a probability 
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space on which we define a ra-dimensional Brownian motion B = ( Bt)o<t<T with 
integer rn > 1. Let us denote by F = {Ft, 0 < t < T} for fixed T > 0, the 
natural filtration generated by process B and augmented by A/p the P-null sets, 
i.e. Ft = cr {B s , s < t} V A/p. 

Let V be the tr-algebra on [0,T] x 1! of .Ft-progressively measurable sets. Let 
p G [l,oo) be real constant and t G [0,T] be fixed, We then define the following 
spaces: C p = {£ : Fj-measurable and R m -valued random variable s.t. E[|£| p ] < oo}; 
Sf T = {^p = (tp s )t<s<T- P-measurable and R m -valued s.t. E[sup sg [ t T ] |vj s | p ] < oo} 
and h(? tT = {y = (<Ps)t<s<T '■ 'P-measurable and R m -valued s.t. E[(/ t T |vj s | 2 ds)5] < 
oo}. Hereafter, T and T are simply denoted by and 
The following assumptions are in force throughout this paper. Let <r be the func¬ 
tion defined as: 

a : [0, T\ x R m — x R mxm 

which satisfies the following assumption. 

Assumption 2.1 (i) a is uniformly Lipschitz w.r.t x. i.e. there exists a constant 

Ci such that, Vi € [0, T ], Vx, x' G R m , |cr(i, x) — <r(i, x')\ < C\\x — x'\. 

(ii) cr is invertible and bounded and its inverse is bounded, i.e., there exits a con¬ 
stant C a such that\/(t,x ) G [0,T] x R m , |<r(i,x)| + |cr _1 (i,a;)| < C a . 

Remark 2.1 Uniform elliptic condition. 

Under Assumptions 2.1, we can verify that, there exists a real constant e > 0 such 
that for any ( t,x ) G [0,T] x R m , 

e.I < a(t,x).a T (t,x) < e -1 ./ (2.1) 

where I is the identity matrix of dimension m. 

Suppose that we have a system whose dynamic is described by a stochastic dif¬ 
ferential equation as follows: for (f, x) G [0,T] x R m , 

( X t,x = x + T a(u, X^)dB u , s G [f, T \; 

l Jt (2.2) 

[x‘- x =x, sG [0,i], 

The solution X = (X*’ x ) s <t exists and is unique under Assumption 2.1. (cf. 
Karatzas and Shreve 1991 [20], p.289). We recall a well-known result associates to 
integrability of the solution. For any fixed (f, x) G [0,T] x R m , p > 2, it holds that, 
P-a.s., 


E 


sup |A'*’T 

.0<s<T 


<c(i + \x\n, 


(2.3) 


where the constant C is only depend on the Lipschitz coefficient and the bound 
of a. In addition, property (2.3) holds true, as well, for the expectation under the 
probability which is equivalent to P. 
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For integer n > 1, we first present the following Borelian function as the terminal 
coefficient of the n-dimensional BSDE that we considering: 

g i :R m —* R, i = 1,2 


which satisfy 

Assumption 2.2 The function g l , i = 1,2, are of polynomial growth with 
respect to x, i.e. there exist constants C g and 7 > 0 such that 

\g i { x )\ < C g (l + |x| 7 ),Vx £ R m , for * = 1,2,..., n. 

Now, we consider Borelian functions Hi,i = 1,2,..., n, from [0, T] xR m xR“xR" m 
into R as follows: 


Hi(t, x, y 1 ,..., y n , z 1 


i = 1 , 2 ,..., n 


which satisfy the following hypothesis: 

Assumption 2.3 (i) For each (t, x, y 1 , ..., y n , z 1 , ..., z n ) £ [0,T] x R m x R" x 

R nm , there exists constants C 2 , Ch and 7 > 0, such that, for each i = 1, 2..., n, 

\H. l {t,x,y 1 ,...,y n ,z 1 ,...,z n )\ < C 2 (l + |a;|)|^| + C h (l + |a;| 7 + |y*|); (2.4) 

(ii) the mapping (y 1 , ..., y n , z 1 , ..., z n ) £ R n xR“ m 1 —^ Hi(t , x, y 1 ,y n , z 1 ,..., z n ) £ 
R is continuous for any fixed ( t,x ) £ [0,T] x R m . 


The BSDE that we concern in this work is the following: 


Yf = g\X°’ x ) + J" T H t (s , X°’*, Y s \ ..., Y s n , Z\, ..., Z n s )ds - £ Z\dB s . (2.5) 


for i = 1,2From Assumptions 2.2 and 2.3, we know that this is a multiple 
dimensional coupled BSDEs system under Markovian framework, with unbounded 
terminal value and the generator of the ith equation is of linear growth on Z l 
component w by w, or in another word, it is of stochastic linear growth. 


3 Existence of solutions for the multiple dimensional coupled 
BSDEs system 

In this section, we will provide an existence result of BSDEs (2.5) when n = 2 as 
an example. Actually, the case for n > 2 can be dealt with in the same way without 
any difficulties. 


3.1 Measure domination 

Before we state our main theorem, let us first recall a result related to measure 
domination. 
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Definition 3.1 ; L q -Domination condition 

Let q £ (l,oo) be fixed. For a given t\ £ [0, T], a family of probability measures 
{v\(s,dx),s £ [ti,T]} defined on R m is said to be L q - dominated by another family 
of probability measures {vq (s,dx),s £ [fi,T]}, if for any 6 £ (0,T — t{\, there exists 
an application <j>\ : \t\ + <5, T] x R m — > R + such that: 

(i) v\(s, dx)ds= 4> s tl (s, x)vq(s, dx)ds on [ti + 5, T]x R m . 

(ii) Vfc > 1, f> s tl (s,x) £ L q ([t\ + S,T] x [— k, k\ m ; oq( s,dx)ds). 

Lemma 3.1 Let xq £ R m , (t,x) £ [0,T] x R m , s £ ( t,T] and p,(t,x;s,dy) the 
law of Xl ,x , i.e., 


VA € B( R m ), n{t,x\s,A) = P (Xl' x £ A). 

Under Assumption 2.1 ona,for any q £ (l,oo), the family of laws {p,(t,x;s,dy),s£ 
[t, T]} is L q -dominated by {/j,(0,Xq; s,dy),s £ [t, T]} for fixed Xq £ R m . 

Proof See [15], Lemma 4.3 and Corollary 4.4, pp. 14-15. □ 

3.2 High dimensional coupled BSDEs system 

Our main result in this section is the following theorem. 


Theorem 3.1 Let x £ R m be fixed. Then under Assumptions 2.1, 2.2 and 2.3, 
there exist two pairs of V-measurable processes (Y l ,Z l ) with values in R 1+m , i = 
1,2, and two deterministic functions U(t,x) which are of polynomial growth, i.e. 
|C(t, a:)| < C(1 + |a;| 7 ) with 7 > 0, i = 1,2 such that, 


P-a.s.,Vf < T,Yf = fi l (t,Xf X ) and Z l is dt-square integrable P-a.s.; 

Y t = 9\X°t' X ) + [ X° ,x , Y S \Y S 2 , Zl, Z 2 )ds - £ Z\dB a ; (3J) 

Y t = g 2 (X?’ x ) + [ T H 2 (s , X 0 *, Y S \Y S 2 , Zl Z 2 )ds - [ T Z 2 dB s . 

Jt Jt 


The result holds true as well for case i = n > 2 following the same way. 


Proof The structure of this proof is as follows. We first use the mollify technique on 
the generator Hi, to construct a sequence of BSDEs with generators which are of 
Lipschitz continuous. Then, we provide uniform estimates of the solutions, as well 
as the convergence property. Finally, we verify that the limits of the sequences are 
exactly the solutions for BSDE (3.1). 

Step 1. Approximation. 

Let f be an element of C°° (R 2+2m , R) with compact support and satisfies: 


R2+2 


ayly 2 ,z\z 2 )dy 1 dy 2 dz 1 dz 2 = l. 
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For (t, x , y 1 , y 2 , z 1 , z 2 ) G [0, T] x R m x R 2+2m , we set, 
x, y 1 , y 2 , z 1 , z 2 ) 

= / n 4 Ri(s,^„(a;),p 1 ,p 2 ,t? 1 ,g 2 )- 

J R2 + 2m 

£ (nfo 1 - p 1 ), n(y 2 - p 2 ), n ^ 1 - q 1 ), n (* 2 - <? 2 )) dp 1 dp 2 dq l dq 2 , 


where the truncation function (p n (x) = ((xj V (—n)) A n)j=i, 2 ,...,m, for x = 
G R m - We next define ip G C'°°(R 2+2m , R) by, 






l, |y 1 | 2 + | 2 / 2 | 2 + k 1 | 2 + k 2 | 2 <i, 

0, IwT + |^| 2 + |*T + |*T>4. 


Then, we define the measurable function sequence (7?i n ) n >i as follows: 
M(t,x,y l ,y 2 ,z l ,z 2 ) G [0,T] x R m x R 2+2m , 


12 12 

tfinfoabS/ 1 ,?/ 2 ,* 1 ,* 2 ) = ip(—, ~, — )Hi n (t,x,y 1 ,y 2 ,z 1 ,z 2 ) 

n n n n 


We have the following properties: 


(a) Hi n is uniformly lipschitz w.r.t (y 1 ,y 2 , z 1 , z 2 )-, 

(b) |-Km(t, Xjy 1 ,y 2 , z 1 , z 2 )\ < C 2 ( 1 + |^„(x)|)|x 1 | + (7^(1 + |pn(x)| 7 + jp 1 !); 

(c) |-ffin(t, 2 / 1 ! 2 / 2 , z 1 , z 2 )| < c„, for any (t,x,y 1 ,y 2 ,z L ,z 2 )\ 

(d) For any (t,x) G [0,T] x R m , and K a compact subset of R 2+2m , 

sup \H ln {t, x, t/ 1 , t/ 2 , z 1 , z 2 ~) - Hi(t, x, j/ 1 , j/ 2 , z 1 , 2 2 )| -a- 0, 

(y 1 ,y 2 ,z 1 ,0 2 )6K 

as n —» oo. 

(3.2) 


The construction of the approximating sequence (-ff^) ra >i is carried out in the 
same way. 

For each n > 1 and (t. x) G [0, T] x R m , since H\ n and H 2 n are uniformly lipschitz 
w.r.t (y 1 , y 2 , z x ,z 2 ), by the result of Pardoux-Peng (see [21]), we know, there exist 
two pairs of processes (Y in ’^ x \ Z in '^) G S 2 T ( R) x U\ T { R m ), i = 1,2, which 
satisfy, for s G [t, T], 
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(3-3) 

Meanwhile, the properties (3.2)-(a),(c) and the result of El-karoui et al. (ref. 
[ 11 ]) yield that, there exist two sequences of deterministic measurable applications 
s ln (resp. g 2n ) : [0, T] x R m — > R and 3 ln (resp. 3 2n ) : [0,T] x R m — >■ R m such that 
for any s G [£, T], 

Y s ln;(t ’ x) = g ln (s, X*’ x ) (resp. F s 2 " ;(t ’ x) = c 2 n (s, X*’ x )) (3.4) 

and 

Z\ n ' (t ’ x ) = 3 ln (s, Xl’ x ) (resp. Z 2n ^ = 3 2n (s,Xl’ x )). 


Besides, we have the following deterministic expression: for * = 1,2, and n > 1, 


S m (t, x) = E g\X^ x ) + / ^(8, Xl’ x )ds , V(f, a:) G [0, T] x R 


(3.5) 


where, 

Fin(s, x) = H in (s,x,<; ln (s,x),<; 2n (s,x), 3 ln (s,x),i l 2n (s,x)). 

Step 2. Uniform integrability of (Y ln ’ < ' t ’ x \ Z 1 " ; ( t,x ))„> 1 for fixed (t,x) G [0,T] x 
R m . 

For each n > 1, let us first consider the following BSDE: 

Y s ln = g\xF) + jT C 2 (l + \y n (Xl' x )\)\Zl n \ 

+ C h (l + \ i p n (X t r ’ x )\^+ \Y r ln \)dr~ [ Zl n dB r . (3.6) 

J S 

For any x G R m and n > 1, the mapping (y 1 ,^ 1 ) G R x R m € 2(1 + 
|<p n (X* ,:c )|)| 2 : 1 | + (7^(1 + |(p n (X* ,x )| 7 + I*/ 1 !) is Lipschitz continuous, therefore, the 
solution ( Y ln ,Z ln ) G <S 2 T (R) x "H 2 r (R m ) exists and is unique. Moreover, it fol¬ 
lows from the result of El-Karoui et al. (see [11]) that, Y ln can be characterized 
through a deterministic measurable function ? ln : [0, T] x R m —s- R, that is, for any 
s G [t, T], 


Y s ln = g ln (s,Xl’ x ). 


(3.7) 


Next let us consider the process 

B: = B a - f l [t ,T](r)C 2 (l + |^(X^)|)sign (Z' n )dr, 0 < s < T, 

Jo 

which is, thanks to Girsanov’s Theorem, a Brownian motion under the proba¬ 
bility P” on (UjJ 7 ) whose density with respect to P is St ■= St(Jq C -2 (1 + 
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|v? n (X*’ E )|)sign(Z s ln )l[ tir ](s)dB s ), where for any z = (z*) i= G R m , sign(z) = 
(l[|«*|#o]]Fy)*=l,...,d and £t(-) is defined by 

£{M) := (exp {M t - (M) t / 2})*<t, (3.8) 


for any (,F t , P)-continuous local martingale M = {M t ) t < T . Then (3.6) becomes 
Y} n = g\X^)+j^ C h (l+\tp n (Xp x )\' y +\Yr n \)dr—2} n dB?, s G [t,T]. ( 3.9) 

Then, by using Ito’s formula to e Cht Y t ln , t <T, we know, 

/■ T /- T 
= e c k T gi(jtij _ / C^e^l + | V „(Ap’ x )p)dr + / e Chr Z 1 r n dB?. 

Jt Jt 

Considering (3.7), we have, 


S ln {t,x) = e~ Ght E" 


s^VP^*) + J C h e GhS { 1 + |^(^> a: )P)^ 


•Ft 


where E n is the expectation under the probability P”. Taking the expectation on 
both sides under the probability P n and taking account of <T ln is deterministic, we 
obtain, 


-In 


(t,x) = 


-c h t 


E n 


s c ^g\X^ x ) + jT C h e ChS { 1 + |v?„pC’ x )P)ds 


Since, g 1 is of polynomial growth, and e c ’'* t < 1 for t £ [0,T], we infer that, 
V(t,x) G [0, T] x R m , 


|? ln (t, *)| < CE n 


sup (1 + |A^P) 

t< s <T 


< C E 


St- sup (l + |X s 4 ’T) 

t<s<T 


where the constant C depends only on T, Ch and C g . It follows from a result of 
Haussmann ([16], p.14, see also [15], Lemma 3.1) that, there exists some po G (1,2) 
and a constant C which is independent on n , such that, E[|£t| po ] < C uniformly. 
As a result of Young’s inequality and estimate (2.3), we have, 


S ln (t,x)\ <C P 0 {E[|£ T r]+E 


sup 

t<s<T 


(l + I A* 


. p o'y 

po-i 


}, 


which yields, 


\g ln (t,x)\ < (7(1 + |x| A ) with A = p 0 7/(Po - 1) > 2. 


(3.10) 
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Next by the comparison theorem of BSDEs and property (3.2)-(b), we actually 
have, for any s £ [t, T], 

Y^’ x ) = ^ n {s,Xl’ x ) < Y s ln = q ln (s,Xl’ x ), 

and by choosing s = t, we get that <; ln (t,x ) < (7(1 + |.r| A ), (t, x) £ [0 ,T\ x R m . 
In a similar way, we can also show that for any (t,x) £ [0, T] x R m , q ln (t,x) > 
—(7(1 + |x| A ). As a conclusion, c ln is of polynomial growth on (f, x) uniformly in n, 
i.e. there exists constants C which is independent of n, and A > 2 such that, 

\, ln (t,x)\ <C(1 + M a ). (3.11) 


Combining (3.11) and (2.3), we deduce that, for any a > 1, i = 1,2, 


E 


sup |y; n; ^ x )| 

t<s<T 


< c. 


(3.12) 


On the other hand, by applying Ito’s formula with (Y m '^ t,x ^) 2 and considering the 
uniform estimate (3.12), we can infer in a regular way that, for any t £ [0,T], 
* = 1 , 2 , 


E 



< C. 


(3.13) 


Step 3. There exists a subsequence of ((Y s ln ’( 0 ’ x \ zl n ’^ 0,x ^)o< a <T)n>i which con¬ 
verges respectively to (Yj, 1 , ^i)o<s<T, solution of the BSDE (3.1). Let us recall the 
expression (3.5) for case i = 1, and apply property (3.2)-(b) combined with the 
uniform estimates (3.12), (3.13), (2.3) and the Young’s inequality to show that, for 

1 < q < 2, 


q )ds 


(3.14) 

Therefore, there exists a subsequence {nj,} (for notation simplification, we still 
denote it by {n}) and a £>([0, T])<S>Z3(R m )-measurable deterministic function Ei(s, y) 
such that: 


E 


\F ln (s,X° s ’ x )\ q ds 


< (7E 


(1 + |^ n (A s 0 > x )|)«|Z s ln -( 0 ' x )|« + (1 + \ip n (X® ,x )\' yq + |Y s ln ^°> x )| 


1 } 


< (7{E 

f \zl n *°’ x) \ 2 da 

+ E 

sup |Y s ln;(0 ’ x) | 9 


Jo 


-0 <S<T 


< (7. 


Fin —> Fi weakly in £ 9 ([0,T] x R m ; /r(0, x; s, dy)ds). 


(3.15) 
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Next we aim to prove that (c ln (i, a;)) n >i is a Cauchy sequence for each (t,x) £ 
[0, T] x R m . Now let (t,x) be fixed, 77 > 0, k, n and m > 1 be integers. From (3.5), 
we have, 


<; ln (t,x)-<; lm (t,x) |=E / F ln (s, X\> x ) - F lm (s, X* a ’ x )ds 


< E 


ft+rj 


+ 


E 


E 


/t+rj 

f 

>t+ri 


| F ln (s,Xl’ x )-F lrn (s,X^)\ds 
(■ Finis, Xi •*) - F lm {s,Xt>*)) .1 {xt ..\< k} ds 


(F ln (8,Xt’*) - F lm (s,Xl ’*)) .l uxt ,* l>k} ds 


where on the right side, noticing (3.14), we obtain, 

rt+r) 


E 


| F ln (s,Xl’ x ) - F lm (s,Xl’ x )\ds 


/ \F ln (8,Xi' x )-F lm (8,Xi' x )\*ds }* <Crf^. 

q 

At the same time, Lemma 3.1 associates with the C i- 1 -domination property implies: 


E 


lt+ri 


(Fi n (s,X s ' ) Fi m (s,X s ’ )) •l{|x| ,a ’|<fe} c ^ s 


' R m Jt+r] 

f f 

/ Ft 771 J 


(F ln (s, rj) - F lm (s, v))-l{\rj\ <fc}£t(*> x; s, dr])ds 


(Fin(s , 77 ) - F lm (s, ??)).1{M <k}0t,x(s, 77 ) 77 ( 0 , z; s, dr])ds 


Since ^>t,x{s,i 7 ) £ £«- 1 ([t + 77 , T] x [— fc, /c] m ; /x(0, x; s, dr/)ds), for fc > 1, it follows 
from (3.15) that for each (t,x) € [0,T] x R m , we have, 


E 


Finally, 


/ t+rj 


(■ Fln(s,Xl’ x ) - F lm (s,Xl’ x )) l{| x t,*|< fe }ds 0 as 77, 772 00 . 


E 


Jt+r] 


( Fln{s,X s ’ ) Fl m (s,X s ’ )) -l{|_Y*’ x |>fc}^ S 


< C{E 


't +17 




}^{E[ / |Fi„(s, X‘’ x ) — Fi m (s, Xl ,x )\ q dsj } 9 

L -7t+I7 J 


q— 1 

< Ck~~. 


Therefore, for each (t,x) £ [0, T] x R m , (<j ln (t, a;)) n > 1 is a Cauthy sequence and 
then there exists a borelian application c; 1 on [0,T] x R m , such that for each 
(t-,x) £ [0,T] x R m , lim n _ ) . 00 <j lri (t, :r) = c 1 (i,a;), which indicates that for t £ [0, T], 
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limn-j.oo Y^ n ^°' x \uj) = ? 1 (i, X°’ x ), P—a.s. Taking account of (3.12) and Lebesgue 
dominated convergence theorem, we obtain the sequence ((YJ ln ’( 0,a ^)o<t<T)n>i con¬ 
verges to Y 1 = (? 1 (t, Xt’ x )) 0 < t < T in £ P ([0,T] x R m ) for any p > 1, that is: 


E 



iyln;( 0 ,®) 


Y?\*dt 


as n —> oo. 


(3.16) 


Next, we will show that for any p > 1, Z ln ’^°’ x ^ = ( 3 lra (t, X°’ x )) 0 < t <T)n>i has a 
limit in %j~( R m ). Besides, (Y ln '^°’ x ^) n >i is convergent in <S|,(R) as well. We now 
focus on the first claim. For n,m > 1 and 0 < t < T, using Ito’s formula with 
(Y t ln —Y t lm ) 2 (we omit the subscript (0,x) for convenience) and considering (3.2)- 
(b), we get, 


I Y t ln -Y t lm \ + / \Zl n -Z] m \ 2 ds 


= 2 j\Y s ln - Y s lm )(H ln (s, X° s ’ x , Y s ln , Y 2n ,Z s ln , Z 2n )- 
- H lm (s, X°’ x , Y s lm , F s 2m , Zl m , Z 2 s m ))ds - 2 j\Y s ln - Y s lm )(Zl n - Z] m )dB s 


<cj |r s ln -i; lm | 


(z* + \zt m m + \x. 


0,^1 


+ \Y s ln \ + |y s lm | + (1 + \X° s ’ x \y ds - 2 / (Y s ln - Y s lm )(zy - Zl m )dB s . 


Since for any x,y,z € R, \xyz\ < ±|:c| Q + \\y\ b + ±\z\ c with i + I + i = 1, then, 
for any e > 0, we have, 


| Y t ln - Y t lm | :2 + J \Z] n - Z] m \ 2 ds 

< C£J\\Z?\ + \Zl m \fds + j J\ 1 + |X?-|)‘* 

+ i / T ^ ~ + I / T(l ^ n| + l^ lm |) 2ds ( 3 - 17 ) 

+ | J\ 1 + \x°’ x \)^ds + I Y s ln Y s lm \ 2 ds } 

r T 

-2J (Y s ln -Y s lm )(Zj n - Zl m )dB s . 

Taking now t = 0 in (3.17), expectation on both sides and the limit w.r.t. n and m, 
we deduce that, 


lim sup E 

n,m—> oo 



(3.18) 


due to (3.13), (2.3) and the convergence of (3.16). As e is arbitrary, then the se¬ 
quence (Z ln ) n >i is convergent in 'Wp to a process Z 1 . 
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Now, returning to inequality (3.17), taking the supremum over [0,T] and using 
BDG’s inequality, we obtain that, 


E 


sup | Y^-Y, 1 
0<t<T 


I z\ n - z\ 


l ds 


r_ 

L 2 


£ 

T 


i] 

4 


' sup \Y t ln — Y t lm \ 2 

+ 4E 

' f T \Z] n - F s lm | 2 ds 

L 0 <t<T 


Jo 


which implies that 


limsupE sup | Y t ln 

n,m—>oo *- 0<t<T 


-^rlm |2 


= o, 


since e is arbitrary and (3.18). Thus, the sequence of (F ln ) n >i converges to F 1 in 
<Sj. which is a continuous process. 

Finally, repeat the procedure for player i = 2, we have also the convergence of 
{Z 2n ) n > i (resp. (F 2n )„>i) in U 2 T (resp. S 2 ) to Z 2 (resp. Y 2 = c 2 (-, 


Step 4- The limit process (F/, Zl)o<t<T,i = 1, 2 is the solution of BSDE (3.1). 
Indeed, we need to show that (for case i = 1): 

F 1 (t,X?’ x )=H 1 (s,X° s ’*,Y s \Y s 2 ,Z 1 s,Z 2 ) dt®dP-a.s. 


For k > 1, we have, 


E 


[ T | H ln (s, X°’ x , Y s ln , Y 2n , Zl n , Z 2n ) - Hi(s, Xg’ x , F/ ,Y 2 , Z \, Z 2 )\ds 

JO 

E [ [ T | H ln (s, X°' x , Y s ln , Y 2n , Z l s n , Z 2n ) 

.Jo 

- Hi(s, X®’ x , F s ln , Y 2n , Zl n , Z 2n )I • 

[ T \H ln (s , X s °>", F s ln , F s 2 ", Z a ln , F 2 ") 

.Jo 


E 


- Hi{s, X° s ' x , Y s ln , Y 2n , Zl n , Z 2n )\ • l { | nl „ |+ | r ^| +|zl „| +|z ^|> fc} d S 


E 


[ T | H^s, X°’ x , Y s ln , Y 2n , Zl n , Z 2n ) - Hiis, X° s ’ x , F s \ F s 2 , Zl,Z 2 )\ds 
.Jo 


: = /r + /« + J 3 «. 


(3.19) 


The sequence /", n > 1 converges to 0. On one hand, for n> 1, the point (3.2)-(b) 
in Step 1 implies that, 


| H ln (s,X"’ x ,Y s ln ,Y s 


'-In \r2n ry\n ry2n 


Z 2n ) - H^a, x°’ x , Y s ln , Y 2n , Zl n , Z 2n ) \ 


■ 1 {|v,i"|+|y/»|+|z»"|+|z?"|<fc} 

< C 2 (l + \X°' x \)k + C h { 1 + |X s °’ x p + fc). 
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On the other hand, considering the point (3.2)-(d), we obtain that, 


I H ln (s, I"’ 1 , Y s ln , Y 2n , Zl n , Z 2n ) - Hi(s, X®’ 00 , Y s ln , Y 2n , Z s 1? \ Z 2n )\ 


■ 1 




< sup | H ln ( s , X"’*, yl,y 2 s , z L s ,z 2 s ) - H^s, X"’ x , y\, y 2 s , z], z 2 s 

\vl\+\y 2 s\+\zl\+\zl\<k 

0, 


as n —> oo. Thanks to Lebesgue’s dominated convergence theorem, the sequence /" 
of (3.19) converges to 0 in 

The sequence /£“ in (3.19) is bounded by k2 ^i)/ q with q £ (1,2). Actually, from 
point (3.2)-(b) and Markov’s inequality, for q £ (1,2), we get, 


I 2 < 


c{e[ / (1 + ix,°->i)+ (i + ix;.»r + |y;“i) s <i»]} s 




x E 


1 


Jo 


{\Y s ^\ + \Y^\ + \Z^\ + \Z^\>k} 


ds 


]}' 


pT pT 2 

<c\ E / \z] n \ 2 ds + E / (1 + 1 x*' x \)*%ds 
^ '-Jo - 1 Jo 


E 


'J \Y s ln \ 2 ds] +E[^ (1 + |X s ^|)^ s ]} + 


{e[/ 0 T \Y s ln \ 2 + |yff + \z^\ 2 + |Z*f da] } 


(fc 2 


C 

— 2(1-1) ' 
k 1 


The last inequality is a straightforward result of the estimates (2.3)(3.12) and (3.13). 


The third sequence / 3 1 , n > 1 in (3.19) also converges to 0, at least for a subse¬ 
quence. Actually, since the sequence (Z ln ) n > 1 converges to Z 1 in then there 
exists a subsequence {Z lnk )k> 1 such that it converges to Z 1 , dt ® dP-a.e., and 
furthermore, sup fc>1 |Z t lrlfc (w)| £ On the other hand, (Y lnk )k> 1 converges to 
Y 1 > 0, dt®dP-a.e.. Thus, by the continuity of function H\(t, x, y 1 , y 2 , z 1 , z 2 ) w.r.t 
(y 1 , y 2 , z 1 1 z 2 '), we obtain that 

H\(t X^' x Y t lnk Y^ nk Z^ nk Z^ nk ) 

^ Hi{t, X°’ x , Y t \Y 2 ,Zl,Z 2 ) dt®dP- a.e. 

In addition, considering that 

sup|F 1 (t,X+,y+,F t 2 " fc ,Z+,Z+)| £ n q T ( R) for 1 < q < 2 , 

k> 1 
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which follows from (3.14). Finally, by the dominated convergence theorem, one can 
get that, 


Hi (t, X°' x ,Y t lnk , Y t 2nk , Zl 


k—too 




Z 2u k ) 

,Y t \Y t 2 ,Zl,Z 2 ) in H q T , 


which yields to the convergence of /J in (3.19) to 0. 

It follows that the sequence X^’ x , Yj n , Yf n , Z^ n , Zf n )o<t<T)n>i con¬ 

verges to (Hi(t, X®’ x , Yj, Yf, Zj, Zf))o< t <T in £ 1 ([0,T] x Q,,dt<g>dP) and then 
Fi(t,Xt’ x ) = Hi(t,Xt’ x ,Yt,Y?,Zl,Zf),dt®dP-a.e. In the same way, we have, 
F 2 (t, X®' x ) = H 2 (t, X^’ x ,Yj,Yf,Zi,Zf), dt ® dP-a.e. Thus, the processes (Y 1 , Z z ), 
i = 1,2 are the solutions of the backward equation (3.1). 
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